QCD perturbation theory at large orders with large 
renormalization scales in the large limit. 
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Abstract. We examine the QCD perturbation series at large orders, for different values of 
the 'large (Sq renormalization scale'. It is found that if we let this scale grow exponentially with 
the order, the divergent series can be turned into an expansion that converges to the Borel 
integral, with a certain cut off. In the case of the first IR renormalon at 2/Pq, corresponding 
to a dimension four operator in the operator product expansion, this qualitatively improves 
the perturbative predictions. Furthermore, our results allow us to establish formulations of 
the principle of minimal sensitivity and the fastest apparent convergence criterion that result 
in a convergent expansion. 
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1 Introduction 



It has been known for many years that the QCD perturbation series suffers from factorial 
divergencies At the same time, we know that the perturbation theory itself does not tell 
us the whole story and one needs to consider non-perturbative contributions. Because of the 
asymptotic freedom, we expect these contributions to become more important for lower ener- 
gies. This is confirmed by the operator product expansion (OPE) [2], the general framework 
that is used to parameterize the non-perturbative contributions in several power corrections, 
corresponding to several condensates 0. In fact, one can show that for one source of diver- 
gencies in the perturbation series, the so called IR renormalons, there exists a one to one 
correspondence with these condensates |^ |3] . This is traced back to the Borel insummability 
of the perturbation series, leaving ambiguities in the definition of its Borel sum, that can be 
exactly compensated by different values of the (non-perturbative) OPE condensates. Another 
crucial problem with the Borel sum, is that the Borel integral is expected to diverge at infinity 
pp. However, for now we will disregard this problem and assume that the integral converges 
at infinity, as is the case for several series in the large Pq limit. This makes it still possible 
to define a Borel sum, by choosing a certain prescription for the Borel integral. At large 
energies, the total result for a physical quantity is then supposed to consist of a Borel sum of 
its perturbation series, combined with some OPE power corrections, originating from certain 
values of the non-perturbative condensates. 

For a large class of observables, the dimension four gluon condensate is the lowest dimen- 
sion condensate in the OPE ^. So the validity of the perturbative expansion for a physical 
quantity, depending on one external momentum Q, is intrinsically limited by a Q~'^ term. 
However, the divergent perturbation series gives a stronger limit, coming from the first UV 
renormalon, positioned at — l//3o in the Borel plane. (Incidentally, for the definition of the 
Borel tranform we use in section 2, the first UV renormalon actually lies at -1.) Dealing with 
the series as an asymptotic series, one concludes that, at the 'best' order of truncation, the 
perturbative prediction can be trusted only up to a term. In i,6l, the renormalization 
scale dependence of the uncertainty was analyzed, in the large /3o limit of QCD. It was 
found that for large orders, the size of this term, can be suppressed by increasing the scale, 
or equivalently, decreasing the coupling constant in the perturbative expansion. 

In this paper we will will extend this analysis, by allowing the renormalization scale to 
vary with the order of truncation. We will show that if one lets k (~ In/i^), grow linearly 
with the order of truncation A^, the divergent perturbation series is turned into an expansion 
that converges to the Borel integral, cut off at N/k = at least if x > 1/2- This allows 
us to formulate an expansion that converges to the Borel sum up to a term oc Q"^, roughly 
compatible with the dimension four gluoncondensate. 

In section!^ we examine the large order behavior for the expansion of a generic QCD ob- 
servable with the typical singularity structure in the Borel plane. We arrive at an asymptotic 
formula (|21|) . describing the large order truncations for general values of x- This formula is 
then applied to different toy expansions in section |21 First we consider the hypothetical situa- 
tion, were one has only UV renormalons, then we include the IR renormalons. Furthermore we 
will use formula (|21j) to establish formulations of the principle of minimal sensitivity (PMS) 
and the fastest apparent convergence criterion (FACC) that lead to a convergent expansion. 

^We consider massless QCD, so the quark condensate does not contribute to the OPE. 
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2 The 'scale' dependence at large orders 



Let us start by writing the general expression of the series expansion for a physical quantity 
TZ, depending on one external scale Q, 

Tl{ao) = ao{ro + riao + r2al + . . .) , (1) 

with Oo = Poas{Q) (~ 1/ ln((5^/A^) for large Q). We will examine how the expansion behaves, 
if we reexpand in some new coupling a{k): 

In the large /?o limit, this variation corresponds to a change in the renormalization scale (RS) 
Q Of course nothing prevents us from considering the same variation for real QCD 

expansions; but it then no longer corresponds to a change of the RS. 

The truncation [7^(ao)]Af+i = ao(''o + '^lOo + • • • + rNCiQ), is conveniently rewritten as a 
Borel integral: 

/ dte «o[F(t)]^, (3) 
Jo 

with 

[F{t)]N=ro + nt + r2- + ... + rNj;^^. (4) 
Under the variation the truncated expansion becomes 

[7^(a(fc))]^+l =a(A;)[ro(A:)+ri(A;)a(A;) + ... + r;v(fe)a(A;)^] = / [e'=*F(t)]^ , (5) 

Jo 

with [e^^F{t)]N the obvious generalization of (0J). One now generally assumes that the series 
fo + ^li + r2t^ /2 + . . . has a finite (nonzero) radius of convergence, therefore it defines a 
function F{t), which is analytical around the origin. In fact, it is established, through the use 
of several techniques, that F[t) has several types of singularities on the real axis ^[71 1^ |3] 
at a nonzero distance of the origin: on the negative axis one finds the UV renormalons at the 
points —1, —2, . . .; on the positive axis there are IR renormalons at the points (1), 2, . . . and 
instanton-anti instanton singularities at the points 4/3o'7r, S/^qti", .... With the number of active 
flavors being 3-6 (/?o = (33 — 2Nf)/12Tr), the singularities near t = are all of the renormalon 
type. In the following we will consider the situation described in the introduction, with no 
IR renormalon at t = 1. 

One can now make use of the Cauchy integral theorem for analytical functions (8] , to write 
[e'^*F(t)]iv as a contourintegral along the contour Cq around the origin, excluding the point t 
and the singularities on the real axis (see Figure 1): 

[e''F{t)y = y^{e''F{t)r'' 



\t=o 

n=0 
N , 

V / dze'^^Fiz) 



/ dze^'F{z) 

JCn 



27rz Jco z — t 

2mJco z-t ^z' ^ ' 
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Figure 1: The contour Cq, equivalent to the contours Ci/^, Cuv and Cir, if the residu at 
t < 1/x is picked up. 



For the last equality, we used the analyticity of e^^ F{z) / {z — t) inside Cq. 

We are interested in the asymptotic behavior of if we put k = Nx, the integrandum 
can be cast in the form (...) x e^^'*^^"''^^-', which calls for the saddlepoint technique jS,. So 
we deform Cq to a circle Ci/^ with radius appropriate for a saddlepoint evaluation. If 
t is lying inside the circle, one has to pick up the residu e^*F{t), coming from the factor 
\/{z — t) in the integrandum. As shown in the figure, for 1/x larger than one (two), one 
also has to add the contour Cuv (Cir), in order to exclude the singularities and branch lines. 
For simplicity, we will limit ourselves in the following to pole singularities (without branch 
lines) at the points Zj = 2, 3 . . . on the positive axis and at the points = —1, —2 ... on the 
negative axis, arriving at 



1 



N 



ei- -t)e'''F{t) 
X 



1 

27rf 



:t 



N+l 



dz e^^x^-i-^^) 
z - t 



F{z) 



1 pN{ziX-\riZi) 
Zi X Zi{Zi L) 



(7) 



1 „-Ar(b'|Y+ln|2'|) 

+ E^(--i^:i)(-i)"t""^ ^ 

X 



t] 



ReSz=AFiz)] , 



with 6{x) the unit stepfunction. 

One can check that the singularities at t = Zi of the term 9{l/x ~ t)e^^F(t), are exactly 
cancelled by the corresponding terms in the summation over the singularities at the positive 
axis. Therefore, when performing the t integration of the separate pieces in 0, one can use 
any prescription, as long as the same prescription is used for every piece. We will use the 
Cauchy principal value prescription, defined as: 



PV 



2 Jo 



+ 



(8) 



4 



Before performing the t integration, let us first work out the z integration in the second 
term of (O. At the saddle point, z = one has a a zero derivative for N{xz — Inz). The 
second derivative is +Nx^, so the path of steepest descent crosses the saddle point orthogonal 
to the real axis. Furthermore for z = e*^/x, 

Re{N{xz - Inz)) = iV(cos 6^ + Inx) <iV(l + lnx), (9) 

thus Ci/^ is indeed an appropriate choice for the saddle point contour. The integral is conve- 
niently written as: 

With the reparametrization, 9 = {2/N)^/'^x, this expression is readily expanded in powers 
of l/N^^"^. However, one should be cautious for t ~ 1/^: in that case, the denominator in 
(fTU)) can not be expanded as 

^ +...). (11) 



l-Xt' \Nl-xt 

Later on we will see that the dominant contribution of the t integration of H1U|) in ©, comes 
exactly from this critical region, therefore we have treated A = (iV/2)i/^(xt - 1), as an order 
1 term in the expansion of ()1U() . arriving at: 

i.,A^+l,A^(l+lnx)^^(i)(;„(A) + y|/i(A) + 0(1)) , (12) 

with 

/+00 p-a;2 
dx— — , (13) 



+ 00 



oo 



2 3 F'(i) , ^2 



/i(^) = / dx —-ir[^-x—^) + , . (14) 



x + iA\3 xF{-)^ (a; + iA)2 2 

X 



One can verify that the odd function /o(A) has a discontinuity at A = {t = l/x)i /o(0±e) = 
ibvr, exactly compensating the discontinuity of the term 6{x^^ ~ t)e^^F{t) in Q. Obviously 
the discontinuities coming from the other step functions 9{x^^ ~ also compensated 

by the integration along Ci/^. In principle one could write 

and incorporate these singularities in the same way as we did for the singularity at z = t. This 
complicates the calculations and only modifies formula (|12j) for x^^ ~ ^i/l^'il + 0{1/N^^'^), 
effectively spreading the discontinuous unit step of 6{x'^ ~ ^i/l^'il) ©> over a region 
X-' = z,M\±0{l/Ny^). 

We are now ready to perform the t integration of ((T)) in (jSJ. First, we will consider the 
integration over the saddle point contribution (|12() . Notice that the integrandum can be cast 
in the form (...) x e^(''^*~^*\ with the opposite power of the exponent as for the z-integration. 
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so again we have a saddle point at but now with the positive axis as the path of steepest 
descent. With the reparametrization 

1 [Yx , , 

the integral reduces to: 

2ttx V N J-(Ar/2)i/2 V A* L V A* 



/•+00 „ , 12 \ / 1 \ 

dxe-^ (/o(x)[x(l - — ) + 3^'] + /i (1 + . 



In the last step we used fo{—x) = —fo{x) and fi{—x) = fi{x). With the introduction of 
some auxiliary integrations, the integral is calculated exactly, arriving at 

Each term in Q, coming from a renormalon or instanton-anti-instanton at Zj (z-), con- 
tributes a term 

(...)x / dte-'^^^''''^- (19) 

Jo Zi-t 

to the t integration ©. With the same reparametrization this is easily expanded (for 
1/x ^i) as: 

After collection of all the terms, we finally arrive at the asymptotic formula for the trun- 
cated expansion in a{k) (k = Nx), of an observable TZ{ao) with Borel transform F{t): 



1 



[7^(a(fc))k+l = Pvf^dte-^F{t) + ^'-^({l--^)F{^)-^)(l + 0{^)) 

1 1 iV(2x-l-ln(2x)) 1 X 

+«(- - 2) e~fle.„,[F(.)| ViV (' + '5(iv') f^^' 

+e(- - 1) (-l)^e-i^fle.,=_i[F(z)] J (1 + 0(-; 

X X(X + 1) V A/ V A/ 

We have only kept the contribution of the singularity closest to the origin on the posi- 
tive/negative axis at the point 2/ — 1; for all values of Xj this singularity dominates the 
contributions of the other singularities on the positive/negative axis by a factor c^, with 
01. 

From the discussion on the validity of H12|) , we know that the asymptotic formula is invalid 
at the 'steps' of the Q functions, x"^ = 2/ - 1 ± ©(l/A^^/^)^ 

Also note that the approximations ()18|) and 1)20(1 . for the t integration are only valid for 
A^ ^ l/(xao)- For the relevant values of x and ao that we will use later on, this means 
A^ ^ 10. Therefore, the formula ()21|) can certainly not be used to examine for instance the 
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NNLO {N = 2) result. It should be treated as an asymptotic formula, valid for (very) large 
N, assessing the convergence/divergence of the truncations in the limit N ^ oo for a certain 
value of X- 

We end this section by quoting the result for the asymptotic formula, in the case of a 
Borel transform with general singularities. One can verify that the renormalon contributions 
coming from the z integrations along Cjr and Cuv change to: 

1 1 sminb) r. 1 gAf(2x-l-ln(2x)) [2^ I 

% - ^> A-liir(-.) N^- - + (22) 

X TT XU + 1) V iV V iV 

for the singular behavior. 



Notice that the exponentials, responsible for the large order behavior, have the same argu- 
ments as in (|21|) . For x > 1/2, the first two terms in (|2H) stay the same, and we have 
essentially the same features for the x dependent large order behavior as in the case of pole 
singularities. 



3 Applications 

3.1 The UV renormalons 

It is instructive to consider first the fictitious case of a series that has no IR renormalons 
at all. In that case the IR renormalon term in the asymptotic formula ()21|1 should be left 
out, leaving us with the UV renormalon term as the only possible source of divergence. For 
X < X*) with X* ^ 0-28 the solution of 

X + l + lnx = 0, (24) 

this UV term diverges exponentially at large orders while it it disappears exponentially for 
X > X*- So, only in the latter case the expansion is convergent: 

[n{aiNx))]N+i ^='^ pdt e-^oF{t) . (25) 
Jo 



This is illustrated in Figure 2, where we plot some approximants [TZ{a{Nx))]N+i as a 
function of at = 1/2 for the Borel transform F{t) = 1/(1 + 1). This specific model was 
examined some time ago by Stevenson in the context of the PMS. One clearly observes the 
behavior predicted by (|21() : for x > X* the expansion converges towards the Borel integral 
cut off at 1/x, for x < X* it diverges. The alternation of the divergence originates from the 
factor (—1)^ in the UV term of the asymptotic formula H21|) . At the transition region x ~ X*, 
one has an extremum for odd orders and an inflection point for even orders. 
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Figure 2: The truncations [Tl{a{Nx))]N+i, for N =3[grey dashed], 9 [dashed], 49[grey] (left) 
and N =4, 10, 50 (right) of the series with Borel transform F{t) = 1/(1 + 1), at oq = 1/2, as 
a function of x- Also plotted: the Borel integral cut off at 1/x [full line]. 



From ()21() one easily shows for a general F{t) (with no IR renormalons) , that for large 
odd orders an extremum occurs at: 

Hn) = Nx ^=°° Nx. + , (26) 

that is if ReSz=-i[F{z)]/F{l/x-k) > 0; in the opposite case the extremum occurs for even 
orders. Furthermore, the substitution of H26() in (PT|) shows that, although x ~^ X*-, the 
expansion with k at the extremum k still converges: 

[n{a{k))]N+i '^=°° [~dte~^oFit), (27) 
Jo 

in agreement with |2] for F{t) = 1/(1 + 1), also confirmed in figure 2. 

With ReSz=-i[F{z)]/F{l/x-k) > (< 0), there is no extremum at A; ~ Nx*, for large even 
(odd) orders. If 

1 ^'(— ) 

->-^ + 2x., (28) 
formula (PT|) instead predicts a zero of the second derivative at: 

3 1 Ar 

HN) = Nx ^=°° Nx. + :r^^ . (29) 

Generally, the condition ()28() is fulfilled for ao small enough; for F{t) = 1/(1 + t), it is met 
for oq < 2.9. Again, we have 

1 

[n{a{k))]N+i '^=°° [""^dte-^Fit), (30) 
Jo 

so the PMS converges to the same value both for odd and even orders, at least if we take kpMS 
at the extremum/inflection point in the neighborhood of Nx-k- This is in disagreement with 
the analysis of Stevenson [S] on F{t) = 1/(1 + 1); he predicted the PMS to be biconvergent, 
with an extra term oc e~^^^^*°'°\ for even orders. 
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Figure 3: The approximants [D{a{Nx))]N+i of the large Pq Adler D-function, for 
=3 [grey dashed], 9 [dashed], 49[grey] (left) and N =4, 10, 50 (right), at ao = 0.3, as a 
function of x- Also plotted: the Borel integral cut off at 1/x [full line] and the Borel sum 
[small full line]. 



In the QCD phenomenology one is generally interested in the Borel sum, i.e. the Borel 
integral cut off at infinity (with a certain prescription in the case of singularities at the positive 
t axis). Ordinary perturbation theory, with a fixed value of k, gives a divergent expansion 
with an error cx e"^/"'' ~ Q~^, at the 'best' order. For oq not too large, we can estimate 

t /"OO i I 

r'*dt e~~^F{t) ^ dt e~~>F{t) + 0{e~~>). (31) 
Jo Jo 

Thus our results imply that, with an appropriate choice of k, in the absence of IR renormalons, 
one can approach the Borel sum much better, up to a correction cx Q~'^l^* ~ ■ One way 
to achieve this is to take, for each order, k equal to A^x*- (It can be checked from H21|) that 
this choice for /c, still gives a convergent expansion.) Another way is to apply the PMS, as 
formulated in this subsection. 



3.2 Including the IR renormalons 

Things change when the IR renormalons come into play. Now, the IR renormalon term in 
(|21|) diverges exponentially for x < 1/2; so the expansion only converges for x > 1/2- We 
have performed explicit calculations on the large /3o limit of the Adler D-function. In this 
limit the Borel transform is calculated exact IvpH]: 

32 1 ^ {-\Yv 

with, as expected, singularities for t = —1, —2, . . . and t = 2,3, . . .. Also note that the Borel 
integral converges at infinity, for oq > 0. In order to keep the calculation time under control, 
we have truncated the p-summation at p = 100. In Figure 3, we plot several approximants 
[D{a{Nx))]N+i, as a function of at ao = 0.3, together with the Borel integral cut off at 1/x 
and the Borel sum. We clearly observe the convergence to the cut Borel integral for x > 1/2 
and the divergence for x < 1/2- At the latter region, there is a turnover for odd orders. This 
occurs when the UV renormalon term starts dominating the IR renormalon term in l|21() at 
X = In 2/3 ~ 0.23 for large orders. 
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Due to the singularity of F£){t) at t = 2, the cut Borel integral itself also diverges at 
X ~ 1/2. (For the sake of clarity, we have only plotted the Borel integral for x > 1/3, thereby 
omitting all the other divergencies at x = 1/3, 1/4, . . ..) Therefore, the best value for Xi in 
the sense that the corresponding expansion approaches the Borel sum as close as possible, no 
longer occurs exactly at the transition between convergence and divergence of the expansion. 
Instead, it is recommended to take x slightly larger than 1/2. If oq is not too large, a good 
choice will be 

X = 1/(2 - cao + 0{al)) , (33) 

where c is some positive constant. Indeed, one easily finds that for a general singularity (with 
b > -1) 

the Borel integral cut at 1/x is approximated as: 

/•l/x /-oo _j_ _J_r , 1 _J_ 

/ dte "oF{t)^PV dte "oF(t)+e a^^ K[C + 0{ao)] + 0{aQ) [l + 0(e "o)], (35) 
Jo Jo ^ J 

with C depending on c: 

C = Pvj dt . (36) 

So we find that, with an appropriate choice for k, one can approach the Borel sum up to an 
error oc e"^/*^" oc Q~^, whereas the ordinary perturbation series, with fixed k, gives a minimal 
error Q~^. In addition, if we restrict ourselves in the OPE to the lowest dimension condensate, 
and to the lowest order term of the Wilsoncoefficient, the error can be compensated by 
a redefinition of the condensate. Indeed, with the singular behavior (jSU, the OPE power 
correction (proportional to the ambiguity of the Borel integral) to the Borel sum, is written 
as: 2 

Ge~^ao''K[l + 0{ao)], (37) 

where G parameterizes the experimentally fitted, non-perturbative condensate. Thus, at 
leading order in gq, the Q"^ error l|35|) is absorbed in the condensate by the replacement 
G ^ G' = G -G. 

From the second column of Tables 1 and 2, one clearly observes that the approximants 
with k = N/{2 — 2ao) indeed converge to the Borel integral cut off at (2 — 2ao), for the Adler 
D-function at oq = 0.2 and oq = 0.4. This gives an error of —0.2% and —8% with respect to 
the Borel sum for N — > oo, whereas the divergent series with k fixed at zero, has a minimal 
error of 5% and 29%. Furthermore, there is also a considerable improvement at finite orders. 

As in the case without the IR renormalons, one can also formulate some condition on the 
approximants [7^(a(A;))]Ar+i, which selects for each order a certain value of k. This results in 
a convergent behavior, if for large orders k — > Nx, with x > 1/2. In this region, a condition 
on the approximants, reduces for large A^, to a condition on the Borel integral cut off at 1/x- 
One can for instance define a PMS criterion as: 



^^[7^(a(fc))]jv+l , 

6k^ lfc=fc/ 

which results in the large order behavior 

"1/Xpa/s 



0, (38) 



[n{a{kpMs))]N+i = dte -oF{t), (39) 

Jo 
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n 

u 




kpMS 


kpACC 


1 


-10.77 


-7.62 


-7.26* 


-7.26* 


2 


5.06 


-1.31 


-3.43 


-2.34 


3 


-5.76 


-1.95 


-1.80* 


-1.80* 


4 


5.39 


-1.26 


-1.10 


-0.93 


5 


-7.38 


-1.12 


-0.75* 


-0.75* 


7 


-17.6 


-0.83 


-0.83* 


-0.41* 


10 


150 


-0.63 


-0.24 


-0.22 


25 


-6x10^° 


-0.35 


-0.048 


-0.0070 


50 


8x10^2 


-0.26 


-0.068 


-0.034 


100 


5x10^1 


-0.21 


-0.088 


-0.055 


oo 


iboo 


-0.16 


-0.12 


-0.095 



Table 1: Deviations in terms of percentage of the approximants [D{a{k))]N+i for the Adler 
D-function at oq = 0.2, with respect to the Borel sum, for different values of k. For the 
N ^ oo limits, we took the cut Borel integral, with the cuts determined according to the 
text. 



N\k 





N/{2 - 2ao) 


kpMS 


kpACC 


1 


-33.5 


-21.7 


-21.6* 


-21.6* 


2 


28.8 


-12.1 


/ 


-15.2 


3 


-56.3 


-12.9 


-17.3 


-9.83* 


4 


119 


-11.8 


-16.r 


-8.73 


5 


-283 


-11.4 


-10.8 


-4.96* 


7 


-2685 


-10.7 


-11.5* 


-1.94* 


10 


1.8x10^ 


-10.0 


-14.1 


-2.30 


25 


-2x10^^ 


-8.89 


/ 


-3.18 


50 


9x10^^^ 


-8.40 


/ 


-4.14 


100 


7x10^21 


-8.14 


/ 


-4.71 


oo 


iboo 


-7.85 


/ 


-5.37 



Table 2: The same as in Table 2, now for ao = 0.4. 
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where 

For small values of a^, with the singular behavior one has xpms = 1/(2 — (l + &)ao + 
O(ao)), leading to a error with respect to the Borel sum. However, Eq. (|4U|) does not 
necessarily have a solution for all values of oq- This is the case for the Borel transform Fu{t) 
of the Adler D-function when oq > 0.34. Evidently, one will then no longer find an inflection 
point for large orders, rendering the PMS criterion 1)38(1 useless. Note that the more common 
PMS criterion, with kpMS defined at an extremum, will generally have no solutions at large 
even/odd ^ orders. Furthermore, taking k at the extremum for large odd/even orders will 
not result in a convergent expansion, since this extremum occurs at x ~ In 2/3 < 1/2, as 
mentioned before. 

We stress that our analysis is only valid for large A^. It may well happen that ()38() has no 
solution at certain finite orders, although xpms (|4()|1 exists, implying a solution for (j^H]) at 
large orders. Therefore, at finite orders, the PMS should be applied with a certain amount 
of flexibility. For instance, in the case of the Adler D-function at ag = 0.2, we find that ()38|) . 
has no solution at the odd orders N < 21. At the orders N = 1,3,5, we have taken kpMS 
at the extremum of [D{a{k))]N^i, indicated by a * in the third column of Table 1. At the 
orders N = 7, . . . , 21, kpMS was taken at the (non-zero) minimum of 

5^[Dia{k))]N+i . 

6k^ I' i*^^ 

indicated by a o. (This minimum does not exist for N < 7.) For large orders, one can observe 
the convergence to the Borel integral cut off at 1/xpms- The improvement at finite orders is 
at least equally large (except for = 2) as for the expansion with k = N/{2 — 2ao). 

For oo = 0.4, we find that no value of k can be identified as kpMS^ when A^ > 16, not even 
if we relax the condition, as we did for the low odd orders at oq = 0.2. This is in agreement 
with our analysis, since XPMS does not exist for ao = 0.4. 

We have also considered a FACC, defining kpACC at a minimum of 

[7^(a(A:))]^+l - [7^(a(A:))]^ _ «(fc)^+V^(A;) 

' [lZ{a{k))]N+i ' ' [n{a{k))]N+i ' ■ ^ ' 

From (|21|) one easily shows that the FACC expansion converges 

[n{a{kFACc))]N+i = dte -oF{t), (43) 

Jo 

with Xfacc at the minimum of 

1 e-V(xao)_p(i/^) 



(44) 



For small values one finds again xfacc = 1/(2 — (1 + &)ao + C'(ao)^), thus at leading order in 
ao, the correction equals the one for the PMS expansion. One can show that (if 6 > — 1 
in (EU) a sufficient condition for XFACc{> 1/2) to exist is: F'{0) < 1/ao. For the Adler 

^depending on the relative sign of the IR and UV terms in H21^ 
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D-function, this condition is fulfilled for any positive value of oq, so in contrast with the 
PMS, the FACC will now have a large order solution, both for oq = 0.2 and oq = 0.4, as can 
be seen from the last column of Tables 1 and 2. At large orders, we see the convergence of 
the FACC expansion to the Borel integral cut at 1/xfacc- At low orders, the improvement 
is generally larger than for the PMS or for k = N/{2 — 2ao). At low odd orders {N < 23, 11 
for oo = 0.2, 0.4) kpACC was taken at a zero of (021 • This is indicated again by a *, because 
it is also an extremum of [D{a{k))]N^i^ which follows from 

^[^^"^'j^J^^^ = -{N + l)a(kf^\M{k) . (45) 

4 Conclusions 

The main result of this paper is formula H21|) . describing the asymptotic behavior of trunca- 
tions of the expansion for a generic QCD observable with the typical singularity structure in 
the Borel plane. We have successfully tested the formula on several toy expansions. It implies 
that, with a 'good' order dependent choice for the 'large (3q renormalization scale' k, one can 
recover a convergent expansion from the divergent perturbation series. In the absence of IR 
renormalons, one can make this expansion converge to a limit, differing from the Borel sum 
by a term oc at large Q. For the realistic case, with IR renormalons starting at t = 2 
in the Borel plane, one can obtain a difference oc Q~^, contrasting the minimal error 
of the divergent ordinary (fixed k) perturbation series. To achieve this one can simply take 
k = N/(2 — coo), or, under certain conditions, one can apply some self-consistent criterion 
like the PMS (|38j) or the FACC (|42j) . The latter criteria will generally converge faster. 

Note that in the case of a first IR renormalon at t = 1, a straightforward generalization of 
(|21|). will only show convergence for x > 1, implying a error, again in accordance with 
the lowest dimension OPE condensate. 

For our results to be relevant on a practical level, the convergence has to set in fast enough. 
Although we generally find, that k = N/{2 — cao),k = kpMS oi k = kpACC are also good 
choices at N = 1,2, we stress that one should not expect miracles at these low orders. In our 
opinion, no procedure to determine the optimal renormalization scale and scheme, will ever 
bypass the simple fact that one is using only two or three of an infinite set of coefficients. 
The main thing one should ask for is that the procedure at least results in a convergent 
expansion, with a finite N ^ oo limit, that is (by definition) recovered to arbitrary precision, 
by calculating enough orders. 

Throughout the paper we have worked with perturbation series for which it was possible 
to define a Borel sum, or equivalently, for which the Borel integral converged at infinity. The 
success of the expansions was evaluated by looking how close they approached the Borel sum 
(-|- some OPE power corrections). As mentioned in the introduction, for real QCD, it is argued 
from general principles that the Borel integral in fact diverges at infinity. This obscures 
the picture of the non-perturbative power corrections, naturally arising from the ambiguities 
of the Borel sum, since the Borel sum itself is divergent for any prescription. However, we 
can still use our results, to motivate heuristically the need for a non-perturbative correction 
oc Q~'^ rather than oc (for the first IR renormalon at 2). Indeed, the N ^ oo limit of the 
expansion can be interpreted as an approximant TZ{x) for a physical quantity TZq, dependent 
on one unphysical parameter x- It is then practice to estimate the error A as proportional to 
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small variations of the unphysical parameter. In om' case, this leads to the estimation 

A{x) oc -i-F(-)e-i/('^°^) oc Q^^/^ , (46) 

for large Q and x > 1/2- So the physical quantity TZq, seems to be approached as close 
as possible, by taking x ~^ 1/2, leaving an intrinsic 'error', or non-perturbative correction, 
oc Q~^. To be precise, one should take x = 1/(2 — coq), otherwise the Q independent term 
F{l/x)/x^ ill (|46j) blows up. For the generic singular behavior (|34j) . we then have (for b > —1): 

7^(c') = 7^(c) + e~^^a,'K Jdt^^ + .... (47) 

With the above-mentioned argument in mind, this naturally leads us to the consideration of 
a non-perturbative power correction, 

A(c) = G(c)Kao''e"^ + ... . (48) 

Note that, at leading order, one can fix the c dependence of the 'condensate' G{c), such that 
the total result 'R{c) + A(c) is independent of c, as it should be. 

It remains to be seen, to what extent our results will change if we consider the renormal- 
ization scale variation for real QCD (with nonzero /?i,/325 • • •)■ 
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